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Approximate Hartree-Fock calculations for Hs, CIL, NII2, II2O, H2O+, BHr, Nri4+, 
and Cil4 have bcou used to determine equilibrium angles and internuclear separations. 
The results are in good agreement with the experimental values when^ these exist. 

Data from these calculations are used to test the usefulness of a partition into one- 
electron type terms. It is found that the partitioned energies reflect thc^ influence of the 
other orbitals in such a way as to render a simple interpretation very difhcult. For example, 
the present results predict an equilibrium angle of 119° for 1120+ which is very unexpected 
if the usual picture of a nonbonding ?>i orbital is accepted. 



1. Introduction 

Perhaps the most widely used scheme for the prediction of molecular shapes is the corre- 
lation diagrams proposed by Walsh [1].^ The simplicity and popularity of the scheme rests on 
the assumption that the total energy of the molecule is simply tlie sum of one-electron orbital 
energies. Substantially, this is equivalent to an extended Hiickel calculation as it is usually 
applied [2]. Hoffmann notes that such a procedure can be valid only if there is a slow variation 
with distance for the difference between the nuclear-nuclear repulsion and the electron-electron 
repulsions. 

Keliance upon such a fortuitous cancellation or semiempirical parameterization of the 
problem is not satisfying in two essential respects. In the first place, as noted by Coulson and 
Neilson [3], there are fundamental theoretical objections to any one-electron model. Secondly, 
the predictions of the scheme must be treated more cautiously when applied to changes in 
molecular geometry attendant upon excitation or ionization. 

Coulson and Neilson proposed that the apparent simplicity of the partition of the total 
energy into one-electron terms be retained, nonetheless, for its obvious pedagogic and pre- 
dictive value. The partition which they devolved upon is a natural one based on the Hartree- 
Fock molecular orbital approach. The total energy in the Hartree-Fock approximation can 
be written, for a closed-shell molecule, as 

E=j:,e,+VM, (1) 

i 



I Figures in brackets indicate the literature references at the end of this paper. 
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where e,=£'^+£'^. Y^ is the total nuclear-nuclear repulsion energy, E^ is the expectation 
value for molecular orbital i of all one-electron potential and kinetic energy operators, and 
et is the ^th eigenvalue of the solution of the self-consistent field equations which determine the 
Hartree-Fock function. 

In general, with the exception of Vjv, these quantities cannot be determined accurate^ 
for a general range of molecular cases. The usual method is to expand the Hartree-Fock 
orbitals as a linear combination of atomic orbit als. Solution of the resultant equations is by 
an iterative self-consistent field (SCF) procedure [4]. Accuracy of the variational solutions is 
predicated upon the number and type of basis functions. The total energy of the system is 
determined by the SCF program precisely by eq (1) or its analog suitably modified for the 
open-shell case. 

Implicit in this approach is the assumption that the geometry of a molecule can be cal- 
culated, to a high degree of accuracy, by a self -consistent field one-electron theory which ignores 
certain types of correlation. This assumption is buttressed by both experimental numerical 
work with diatomic molecules [5] and theoretical arguments [6]. 

However, calculations of polyatomic molecules have not been quite as successful. Much 
of this difficulty probably can be attributed to the approximate integrals that were used in 
evaluating the necessary matrix elements. The prediction of the equilibrium angle for AH2 
molecules with a Hartree-Fock calculation, until recently, was quite problematical. 

Higuchi [7] found a rather slow variation of the energy of NH2 with angle in the neighbor- 
hood of the experimental value and a minimum of 105°, but the integrals were not evaluated 
in a consistent fashion. The equilibrium angle calculated by Ellison and Shull [8] was con- 
siderably greater than the experimental value. Using approximate integrals, Krauss and 
Padgett [9] determined the molecular geometry of CH2 which was in poor agreement with 
experimental results. An accurate configuration interaction calculation by Foster and Boys 
[10] did not find any better agreement with the experimental equilibrium configuration for 
either the singlet or triplet ground states of CH2. The total energy, however, is in the neighbor- 
hood of the expected Hartree-Fock limit. This emphasizes the relative unimportance of corre- 
lation in determining the geometry of a molecule. 

SCF results, considerably more accurate than those cited above, have recently been 
published by Moccia [11]. Using a large number of functions in a one-center expansion, he 
obtained equilibrium configurations that are in good agreement with the experimental value 
confirming the impression that the correlation of electrons plays a small role in determining 
the geometry. 

Approximate Hartree-Fock orbitals and energies of the same order as those obtained by 
Moccia can be obtained readily with Gaussian basis sets. Preliminary results for CH4 and 
NH4" [12] were quite encouraging. Qualitatively interesting results have been reported b}^ 
Harrison [13], but the number of basis functions that were used for H2O were insufficient to 
yield a quantitative result. Extending the basis slightly, however, produced both a large 
improvement in the total energy and the calculated equilibrium angle. These results indicate 
that moderately accurate Gaussian SCF results can predict the structures of radicals and ions. 
Calculations for H3, CH2, NH2, H2O, H2O+, BH^, NH4', and CH4 will be reported in this paper. 

The SCF results also yield enough data to test the contentions of Coulson and Neilson. 
It will be seen that the simple one-electron picture is unsatisfactory since the partitioned 
energies reflect the influence of the other orbitals in a complex fashion. Such conclusions 
must be qualified hy observing that accurate Hartree-Fock values have not been obtained 
but these results provide a consistent test with functions that do predict a moderately accurate 
geometry. 

2. Results of Gaussian Calculations 

The SCF calculations follow a standard procedure that has recently been review^ed by 
Nesbet [4]. Since the open shell is of different symmetry from all other occupied orbitals, 
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the Nesbet and Roothaan [14] procedures are equivalent and the SCF equations may be solved 
by the same procedures used for a closed-shell system with a suitably modified A matrix. 
For CH2, ^2~, which can be represented 

the eigenvalues are given by the following expressions 

eu^ = En + y^i \2Jaa^ i r — Ka^^ tt1 + 2j^ „ .^ n — i^icr^,n+c/n^n_— ifn^.n_, 

a. 

where Ki is the expectation value for the ith orbital overall one-electron, potential and kinetic 
energy, operators, Jij is the Coulomb integral between the % and j orbitals, Ki^ is the exchange 
integral between these orbitals, and a=\ or 2. The total energy for CHo, ^2^? i^ then 

with 

Different expressions are required for H3, which is represented by the configuration 

(l(^.)'(l^.). 
The eigenvalues are 

with 

E= eiag+ huu+Eiag+ hEu^, 

A similar result is obtained for NH2, ^Bi, which is represented by the configuration 



where 



and 



€l&2^^'^l''2 ' *^1^2l^2 ' -^ [^^ aa\'\.h2 ^aaj 162J~T~*^1&1 1&2 2^1fti&2; 
a 

a 



The use of Gaussian basis functions is recommended for polyatomic calculations because 
of the ease with which the necessary matrix elements are evaluated (15). However, so many 
functions are required in order to achieve an accurate Hartree-Fock function that in only 
five cases, CH4, NH^, BH^, CH2f 2^), and H3, was an energy obtained that is sufficiently ac- 
curate to warrant an analysis of the energy. For all other cases only the angular variation 
of the approximate energy was considered. 

The parameters for the basis orbitals are given in table 1. The tetrahedral orbital param- 
eters were either scaled from the best H atom results of Longstaff and Singer [16] or were 
varied in a limited manner to minimize the energy. For H2O the choice of parameters was 
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dictated by the desire to compare the present results with those reported by Harrison. The 
parameters for CH2 and NH2 were chosen by reference to various atom results that have 
been obtained by Huzinaga [17]. 

In table 2, the eigenvalues, functions and energies for all of the relevant cases are presented. 
Note that all results are for normalized basis functions. The equilibrium angles and inter- 
nuclear separations were obtained simply from a quadratic fit and only the three solutions in 
the neighborhood of the equilibrium are given. Calculated equilibrium angles and distances 
are compared in table 3 with the experimental values. For the case of BH^, CH2(^2^), and 
H3 a summary of the considerations relevant to determining a bound on the correlation energy 
are presented in table 4. Finally we have tabulated the quantities, et, for CH2, NH2, and 
H2O in table 5. 



Table la. Parameters for Td basis orhitals Nf 
(x,y,z)exp(— ar2) 



Summetry 


Or- 
bital 


Kx,y,z) 


Atom 


type 


H 


B 


C 


N 


fli 


1 
2 
3 

4 

5 
6 
7 
8 
9 
10 
11 
12 

1 
2 
3 
4 
5 
6 
7 


1 
1 
1 
1 
1 
1 

r2 
r2 
r2 

1 

1 

1 

ix,y,2) 
(x,y,z) 
ix,y,z) 
ix,y,z) 

1 
1 




834. 54 
128. 83 
33. 769 
13. 873 
5.697 
2. 0393 
6.08 
0.832 
.307 

9.6 
1.92 

0.448 

.128 


1227. 4 
189. 48 
49. 667 
20. 404 
8.379 
2.999 
9.5 
1.3 
0.48 

15.0 

3.0 

0.7 

.2 


1695. 9 






261. 79 






68. 623 
28. 191 






11.577 






4.144 






13.68 






1.87 






0.69 


H(hi-{-h-^h 
ti 


5.74 

0.775 

.159 


21.6 






4.32 
1. 






0.29 


[H(k-\-h2-h 

-hdMih 

-h2-h 
-h2-\-h 

-h)l 


5.74 
0.775 
.159 





Table lb. Parameters for d^ basis orbitals Nf 
(x,y,z)exp(— ar2) 



Summetry 
type 



Or- 
bital 



H{hi-\-h2) 
bi 

bi 



Kx,y,z) 



^4. 24, 5. 74 

0. 658, 0. 775 
0. 148, 0. 159 



'4.24, 5.74 
0. 658, 0. 775 
0. 148, 0. 159 



1300 
575 
146 
37.3 
9.42 
2.39 
0.608 
.154 
4.78 
0.923 
.178 



4.78 
0.923 

.178 



4.78 
0.923 

.178 



2500 
637.5 
203 
63.5 
20.2 
6.35 
0.7 
.21 
6.89 
1.33 
0.257 



1.; 

0.257 



6.89 
1.33 
0.257 



2500 
1090 
276 
70.4 
17.8 
4.51 
1.15 
0.291 
9.38 
1.81 
0.35 



9.38 
1.81 
0.35 



9.38 
1.81 
0.35 



a The first column of H parameters were used only for the H2O 
calculations. 



Table Ic. Parameters for Dooj^ basis orbitals 

Nf(x, y, z) ea:p(— ar2)(CH2) 











Atom 


Symmetry type 


Orbital 


fix, y, z) 












H 


C 


Ug 


1 






4232. 61 




2 






634. 882 




3 






146. 097 




4 






42. 4974 




5 






14. 1892 




6 






5. 14773 




7 






1. 96655 




8 






0. 496241 




9 






.153313 


l/2(ft,H-ft2) 


10 




5.74 






11 




0.775 






12 




.159 




lo-u 


1 
2 


z 
z 




18. 1557 




3. 9864 




3 


z 




1. 14293 




4 


z 




0. 35945 




5 


z 




. 114595 


moii-hi) 


6 


1 


5.74 






7 


1 


0.775 






8 


1 


.159 






1 


X 




18. 1557 




2 


X 




3. 9864 




3 


X 




1. 14293 




4 


X 




0. 35945 




5 


X 




. 114595 



Table Id. 



Parameters for D 00 /^ basis orbitals Nf (x* 
y, z)ea:p(-ar2)(H3) 



Symmetry type 


Orbital 


/U, y, z) 


H 


hiVA — Vs) - - 


1 

2 
3 
4 
5 
6 
7 
8 
9 
10 
11 

1 
2 
3 
4 
5 
6 
7 


z 

z 
z 


0.8 


UhA+hB) 


22. 683 


he 


3. 6674 

0.972 

.3094 

.1209 

22. 683 


P(, _ __ . _ __ . 


3. 6674 
0.972 
.3094 
.1209 

.8 


if-DX +-Pb") 


.8 


lihA—hB) 


22. 683 




3. 6674 

0.972 

.3094 

.1209 
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Table 2a. BH4~ Eigenvalues^ eigenvector coefficients, and total energy 



Internuclear separation 
(a.u.)* 


2.1 


2.3 


2.4 


Total energy (a.u.) 


-26.9155 


-26.9449 


-26.9451 




lai 


2ai 


U2 


lai 


2ai 


lt2 


lai 


2ai 


1/2 


Eigenvalue. 


-7. 2087 

.0057 
.0398 
.1256 
.2083 
.4266 
.3383 
.0017 
.0164 
-. 0027 
.0004 
.0031 
.0026 


-0. 4625 

-. 0009 

-. 0071 

-. 0198 

-. 0576 

-. 1020 

.0436 

-. 0349 

. 2252 

. 1565 

.0910 

.4946 

.6281 


-0. 1789 

.0118 
.0848 
. 2701 
.1569 
.0980 
.6087 
.9952 


-7. 2676 

.0057 
.0398 
.1256 
. 2088 
.4270 
.3373 
. 0022 
.0157 
-. 0021 
.0002 
.0023 
.0021 


-0. 4434 

-. 0009 

-. 0067 

-. 0183 

-.0579 

-.0966 

.0561 

-. 0383 

.2273 

.1863 

.0847 

.4715 

.6586 


-0. 1765 

. 0109 
.0784 
.2581 
.1581 
. 0942 
.5887 
1. 0300 


-7. 2974 

.0057 
.0398 
. 1257 
. 2089 
.4271 
.3369 
.0023 
.0154 
-.0019 
.0002 
.0020 
.0019 


-0. 4347 

-. 0009 

-. 0066 

-. 0177 

-. 0583 

-. 0943 

. 0618 

-. 0401 

.2284 

. 2007 

.0816 

.4600 

.6740 


-0. 1755 


Eigenvector: 


.0106 


2 


. 0755 


3 


. 2523 


4 _ 


.1597 


5 


. 0928 


6 


.5783 


7 


1. 0459 


g 




9 








10 








11 






12 











a Energy in units of e/2ao= 27.21 eV, 
Length in units of ao= 0.52917X10-! 



Table 2b. CH4 Eigenvalues, eigenvector coeffi,cients, and total energy 



Internuclear separation 
(a.u.) 


1.9 


2.0665 


2.2 


Total energy (a.u.) 


-40.1495 


-40.1668 


-40.1521 




lai 


2ai 


1/2 


lai 


2ai 


1/2 


lai 


2ai 


1/2 


Eigenvalue 


-11. 1691 

.0056 
.0393 
. 1247 
. 2093 
. 4299 
.3341 
. 0023 
.0151 
-.0005 
.0002 
. 0022 
.0004 


-0. 9803 

-. 0010 

-. 0077 

-. 0206 

-. 0729 

-. 1038 

. 0922 

-. 0478 

.2871 

.3000 

.0942 

.4549 

.3612 


-0. 5650 

.0142 
.1005 
.3378 
.2612 
.1021 
. 6327 
.5911 


-1L2121 

. 0056 
.0393 
.1248 
. 2096 
. 4302 
.3334 
. 0026 
.0146 
-. 0003 
.0002 
.0017 
.0004 


-0. 9421 

-.0010 

-.0075 

-. 0195 

-. 0753 

-. 0988 

. 1026 

-. 0526 

. 2839 

.3223 

.0843 

.4322 

.3964 


-0. 5411 

.0136 
.0966 
. 3252 
. 2723 
.0966 
.6148 
.6319 


-11. 2458 

.0056 
.0393 
. 1248 
.2098 
.4303 
.3330 
.0028 
.0142 
-.0002 
.0002 
.0013 
.0004 


-0. 9147 

-. 0010 

-. 0075 

-. 0189 

-. 0773 

-. 0958 

.1095 

-. 0561 

. 2829 

.3381 

.0770 

.4154 

. 4243 


-0. 5234 


Eigenvector: 


.0131 


2 _ -_ 


.0940 


3 


.3158 


4._ - 


. 2803 


5 


.0927 


6 


.5980 


7 


.6670 


8 




9 








10 






11- 








12 

















Table 2c. NH4+ Eigenvalues, eigenvector coefficients, and total energy 



Internuclear separation 
(a.u.) 


1.8 


1.9464 


2.2 


Total energy (a.u.) 


-56.4860 


-56.5038 


-56.4383 




lai 


2ai 


1/2 


lai 


2ai 


1/2 


lai 


2ai 


1/2 


Eigenvalue _ _ _ _ 


-15.9302 

.0056 
.0389 
. 1242 
.2094 
. 4327 
.3317 
. 0021 
.0140 
.0007 
.0002 
.0011 
-. 0001 


-1. 6026 

-.0011 

-. 0083 

-. 0212 

-.0882 

-. 1005 

. 1241 

-. 0605 

.3177 

.3963 

.0874 

.3980 

.2849 


-1. 0425 

.0161 
.1189 
.3815 
.3439 
.1045 
.5706 
.3198 


-15.9444 

.0056 
.0389 
. 1243 
.2096 
. 4329 
.3313 
.0023 
.0136 
.0008 
.0002 
.0009 
-. OOOOi 


-1. 5418 

-. 0011 

-. 0082 

-. 0205 

-. 0917 

-. 0961 

.1330 

-. 0656 

.3158 

.4147 

.0756 

.3810 

.2636 


-0. 9930 

.0155 
.1171 
.3680 
. 3621 
.0968 
. 5528 
.3637 


-15. 9578 

.0056 
.0389 
.1243 
.2098 
.4330 
. 3308 
. 0025 
. 0132 
. 0008 
. 0002 
.0006 
.0001 


-L4487 

-. 0011 

-. 0083 

-. 0198 

-. 0963 

-. 0922 

.1445 

-. 0719 

.3177 

.4409 

.0606 

.3501 

.3113 


-0. 9146 


Eigenvector: 

1 


.0149 


2 


.1145 


3 


.3483 


4 


.3844 


5 


.0848 


6 


.5214 


7 


.4488 


8 




9 








10 








11 






12 
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Table 2d. CH2, 


'Ai, Eigenvalues, 


eigenvector coefficients^ and total energy 






Internuclear 
separation 


100° 


105° 


110° 


i2CH = 

2.1 (a.u.) 

Total energy 

(a.u.) 


-38.8017 


-38.8026 


-38.8020 




lai 


2ai 


3ai 


I62 


lai 


2ai 


Sax 


I62 


lai 


2ai 


Sai 


I62 


Eigenvalue _ _ 

Eigenvector: 

1 


-11.3144 

.0036 
.0062 
.0477 
.2098 
.5740 
.3339 

-.0308 
.0209 
.0027 

-.0003 
.0040 

-.0006 
.0047 

-.0137 


-0.9145 

-.0006 

-.0011 

-.0088 

-.0398 

-.1394 

-.1577 

.4535 

.4495 

.0192 

.1183 

.1232 

.0536 

.3453 

.0560 


-0.3934 

.0003 
.0005 
.0047 
.0209 
.0766 
.0817 
-.2430 
-.4530 
.0739 
.3578 
.4993 
.0331 
.2279 
.2391 


-0. 5623 

.0637 
.3324 
.3119 
.0663 
.5181 
.4027 


-11.3120 

.0036 
.0062 
.0477 
.2098 
.5739 
.3339 

-.0308 
.0209 
.0027 

-.0002 
.0038 

-.0006 
.0049 

-.0137 


-0.9095 

-.0006 

-.0011 

-.0088 

-.0400 

-.1401 

-.1586 

.4561 

.4505 

.0183 

.1137 

.1176 

.0531 

.3439 

.0600 


-0.3859 

.0003 
.0005 
.0045 
.0202 
.0738 
.0789 
-.2350 
-.4361 
.0749 
.3621 
.5157 
. 0323 
. 2209 
.2314 


-0. 5712 

.0635 
.3319 
.3092 
.0667 
.5190 
.3860 


-11.3095 

.0036 
.0062 
.0477 
.2098 
.5739 
.3339 

-. 0307 
. 0208 
.0027 

-.0002 
.0035 

-.0006 
.0050 

-. 0136 


-0.9047 

-.0006 

-.0011 

-.0089 

-.0402 

-. 1408 

-.1595 

.4586 

.4502 

.0175 

.1089 

.1115 

.0525 

.3432 

.0652 


-0.3784 

.0003 
.0005 
.0044 
.0194 
.0710 
.0759 
-.2265 
-.4196 
.0759 
.3663 
.5316 
.0313 
.2138 
.2245 


-0. 5793 
.0634 


2 

3 


.3317 

.3081 


4 

5 

6 

7 


.0672 
.5187 
.3699 


8 








9 








10 








11 








12 








13 








14 















Table 2e. CH2, ^S^, Eigenvalues, eigenvector coefficients, and total energy 



Internuclear 
separation 
(a.u.) 


1.8 


1.95 


2.1 


Total energy 
(a.u.) 


-38.8731 


-38.8936 


-38.8915 




la, 


2.. 


lo-„ 


- 


l<r. 


2<T, 


Uu 


IT 


U, 


2a, 


l<r„ 


TT 


Eigen value.. 

Eigenvector: 

1 


-11.2185 

.0012 
.0093 
.0454 
.1544 
.3583 
.4365 
.1478 
.0032 
.0006 
.0003 
.0C14 
-. 0005 


-0.8915 

-. 0002 

-. 0018 

-. 0087 

-. 0320 

-. 0805 

-.1515 

-. 0562 

.4975 

.2893 

.0590 

.3746 

.1620 


-0. 6886 

.0117 
.0652 
.2322 
.3368 
.0592 
.0835 
.4379 
.3200 


-0. 3973 

.0136 
.0841 
.2755 
.5001 
.3617 


-11.2326 

.0012 
.0093 
.0454 
.1545 
.3584 
.4368 
.1473 
.0031 
.0005 
.0002 
.0011 
-.0004 


-0.8635 

-. 0002 

-. 0018 

-. 0086 

-. 0312 

-.0795 

-. 1472 

-. 0582 

.4974 

.3127 

.0534 

.3563 

.1735 


-0. 6520 

.0109 
.0631 
.2145 
.3472 
.0628 
.0771 
.4261 
.3560 


-0. 3944 

.0136 
.0845 
.2751 
.4954 
.3676 


-11.2456 

.0012 
.0093 
.0454 
.1545 
.3584 
.4371 
.1469 
.0030 
.0005 
.0001 
.0009 
-.0003 


-0.8371 

-. 0002 

-.0017 

-.0085 

-. 0307 

-. 0786 

-. 1447 

-.0592 

.4947 

.3368 

.0487 

.3369 

.1852 


-0. 6183 

.0102 
.0609 
.1998 
.3501 
.0697 
.0718 
.4177 
.3899 


-0. 3920 
.0136 


2 


.0849 


3 


.2759 


4 


.4910 


5 


.3721 


6 




7 








8 








9 








10 














11 














12 


































Table 2f. 


NH2, 


^Bi, Eigenvalues, eigenvector 


coefficients, and total 


energy 






Inter- 




























nuclear 






95° 






105° 










115° 






separation 




























ii;NH = 




























1.95(a.u.) 




























Total 






-55.4496 






-55.4537 










-55.4523 






energy 




























(a.u.) 






























lai 


2ai 


3ai 


I62 


Ibi 


lai 


2ai 


3ai 


I62 


l&i 


lai 


2ai 


3ai 


I62 


I61 


Eigen- 
value 


-15.6019 


-L1341 


-0. 4892 


-0. 6149 


-0. 5049 


-15.6953 


-1.1230 


-0.4709 


-0.6332 


-0. 5005 


-15.5879 


-L1128 


-0. 4531 


-0.6486 


-0. 4957 


Eigen- 
































vector: 
































1 


.0029 


-.0005 


.0003 


.0694 


.0965 


.0029 


-.0006 


.0002 


.0687 


.0962 


.0029 


-.0006 


.0002 


.0682 


.0959 


2 


.0096 


-.0018 


.0008 


.3417 


.4640 


.0096 


-.0019 


.0008 


.3387 


.4622 


.0096 


-. 0019 


.0007 


.3368 


.4603 


3 


. 0385 


-. 0077 


.0036 


.3854 


.6611 


.0385 


-.0077 


. 0033 


.3783 


.6628 


.0385 


-.0078 


. 0031 


.3762 


.6647 


4 


.1452 


-. 0295 


.0133 


.0674 




. 1452 


-. 0297 


.0123 


.0688 




.1452 


-.0299 


.0113 


.0699 




5 


.3752 


-.0892 


.0423 


.4941 




.3752 


-.0899 


.0392 


.4968 




.3752 


-.0904 


.0360 


,4941 




6 


.5361 


-. 1829 


.0841 


.3237 




.5360 


-.1843 


.0780 


.2938 




..5360 


-. 1855 


.0717 


.2686 




7 


.0737 


.4764 


-. 2295 






.0740 


.4799 


-.2140 






.0742 


.4830 


-. 1969 






8 


— . 0416 


. 4438 


-. 3721 






-. 0422 


.4479 


-. 3473 






-.0424 


.4488 


-. 3240 






9 


.0018 


.0164 


.0777 






.0018 


.0150 


.0799 






.0018 


.0135 


.0819 






10 


.0025 


.0966 


.3626 






.0023 


.0888 


.3717 






.0021 


.0808 


.3805 






11 


-. 0067 


.1126 


. 5052 






-. 0062 


.1035 


.5331 






-. 0054 


. 0930 


.5593 






12 


.0011 


.0507 


.0413 






.0011 


.0496 


.0393 






.0012 


.0486 


.0367 






13 


-. 0027 


.2790 


.2625 






-.0032 


.2756 


.2428 






-. 0039 


.2756 


.2245 






14 


.0207 


.0223 


.2239 






.0215 


.0249 


.2169 






.0221 


.0294 


.2100 
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Table 2g. H2O, ^Ai, Eigenvalues, eigenvector coefficients, and total energy 



Inter- 
nuclear 
separation 


96° 


105° 


110° 


E0R = 

1.8 (a.u.) 
Total 
energy 
(a.u.) 


-75.8384 


-75.8453 


-75.8450 




lai 


2ai 


3ai 


I62 


I61 


lai 


2ai 


3ai 


11)2 


l&i 


lai 


2a: 


3oi 


I62 


1^1 


Eigen- 
value 

Eigen- 
vector: 

1 

2 

3 

4 

5 

6 

7 


-20.5836 

.0033 

.0059 

. 0451 

.2016 

.5711 

.3461 

-.0303 

.0184 

.0021 

-.OOOOi 

.0026 

-. 0005 

-. 0010 

-.0068 


-1.3816 

-. 0007 

-, 0011 

-. 0095 

-. 0433 

-. 1601 

-. 1675 

.4886 

.5440 

.0145 

.0820 

.1050 

.0692 

.2099 

-. 0138 


-0. 5794 

.0003 
.0004 
.0038 
.0169 
.0654 
.0626 
-.2007 
-.3502 
.0831 
.3761 
.5205 
.0627 
.2769 
.1693 


-0. 6879 

.0755 
. 3555 
.4205 
.0954 
.4760 
.2426 


-0. 5048 

.1002 
. 4529 
.6696 


-20. 5757 

.0033 

. 0059 

.0451 

.2016 

.5711 

.3461 

-.0303 

.0185 

.0020 

. OOOOi 

.0023 

-. 0005 

-. 0006 

-. 0072 


-1.3097 

-.0007 

-.0011 

-.0095 

-.0436 

-.1611 

-.1685 

.4912 

.5476 

.0132 

.0752 

.0957 

.0676 

.2091 

-.0131 


-0. 5591 

.0003 
.0004 
.0035 
.0156 
.0604 
.0579 
-. 1862 
-.3241 
.0854 
.3856 
.5472 
.0593 
.2491 
. 1668 


-0. 7061 

.0745 
. 3517 
.4122 
. 0982 
.4759 
.2189 


-0.4996 

. 0999 
. 4515 
.6710 


-20.5709 

.0033 

. 0059 

.0451 

. 2016 

.5711 

. 3461 

-. 0303 

.0186 

. 0019 

. OOOO2 

.0021 

-.0006 

-. 0004 

-. 0073 


-1.3638 

-.0007 

-.0011 

-.0096 

-. 0437 

-.1614 

-. 1689 

. 4923 

. 5487 

. 0125 

.0716 

.0905 

.0669 

.2100 

.0127 


-0. 5490 

. 0002 
. 0004 
.0034 
.0149 
. 0578 
. 0554 
-. 1785 
-.3116 
. 0865 
. 3901 
.5594 
.0573 
.2371 
.1650 


-0. 7141 

.0741 
. 3501 
. 4097 
.0993 
. 4740 
.2092 


-0. 4967 

.0998 
.4507 
.6718 


8 












9 












10„ . 












11 












12 












13 












14 

























Table 2h. H20"'", ^Bi : Eigenvalues, eigenvector coefficients, and total energy 



Inter- 






























nuclear 






105° 










115° 








125° 






separation 






























i2oH=1.9 






























(a.u.) 






























Total 






























energy 






























(a.u.) 






-75.4427 










-75.4465 








-75.4462 








lai 


2ai 


3ai 


I62 


I61 


la, 


2ai 


3ai 


I62 


Ibi 


lai 


2ai 


3ai 


I62 


I61 


Eigenvalue 


-21.1628 


-1.8434 


-1.0539 


-1. 1732 


-1.1154 


-21. 1536 


-1.8323 


-1.0372 


-1.1859 


-1.1096 


-21. 1437 


-1.8214 


-1.0211 


-1.1964 


-1.1032 


Eigenvec- 
































tor: 
































1 


.0033 


-.0007 


. 0002 


.0841 


.1080 


.0033 


-.0007 


. 0002 


. 0832 


.1077 


.0033 


-. 0007 


.0002 


. 0825 


.1074 


2 


.0059 


-.0012 


.0004 


.4009 


.5135 


. 0059 


-.0012 


.0003 


.3964 


.5119 


. 0059 


-.0012 


.0003 


.3925 


.5102 


3 


.0451 


-.0100 


.0033 


. 4721 


.6087 


.0451 


-.0100 


.0029 


.4720 


.6103 


.0451 


-.0100 


.0026 


.4754 


.6122 


4 


.2016 


-.0456 


.0149 


. 0607 




.2016 


-.0458 


.0134 


. 0624 




.2016 


-.0460 


.0119 


.0639 




5 


.5713 


-. 1686 


. 0563 


.4126 




.5713 


-.1693 


.0508 


.4067 




.5713 


-. 1699 


.0449 


.3965 




6 


.3460 


-.1781 


.0594 


.0951 




. 3460 


-.1790 


.0537 


.0850 




.3461 


-.1799 


.0475 


. 0762 




7 


-.0303 

.0177 

.0019 

-.0000 4 

.0022 

-.0007 

-.0002 

-.0070 


.5395 
.5581 
.0139 
.0775 
.1084 
.0409 
.1890 
-.0438 


-.2030 
-.2063 
. 0939 
.4339 
.5537 
.0382 
.2065 
.0909 






-.0303 
.0178 
.0018 

-.0000 2 
.0019 

-.0007 
.0002 

-.0073 


. 5424 
. 5597 
. 0125 
.0701 
.0962 
.0398 
.1912 
-.0437 


-. 1837 
-. 1866 
.0960 
.4435 
.5699 
.3500 
.1867 
.0843 






-.0303 
.0179 
.0017 
.OOOOi 
.0016 

-.0008 
.0005 

-.0075 


.5449 
.5604 
. 0109 
.0620 
.0832 
.0388 
.1947 
.0430 


-. 1625 
-.1665 
.0980 
.4527 
.5849 
.0312 
.1667 
.0769 






g 














9 














10 














11 














12 














13 














14 





























Table 2i. 


H3, 'Zt 


: Eigenvalues, eigenvector coefficients 


, and total energy 






Internuclear separation 
(a.u.) 


1.5 


1.6 


1.7 


1.8 


1.9 


Total energy (a. u.) 


-1.58123 


-1.58975 


-1.59296 


-1.59215 


-1.58829 




Ug 


Uu 


l<rg 


l<rn 


l<r, 


l<ru 


la. 


Uu 


u. 


Uu 


Eigenvalue _ 


-.70835 

-.04341 
.01194 
. 06917 
. 24285 
. 32236 
. 12731 
.00861 
. 05011 
. 17421 
. 24290 
.11472 


-0.36524 

. 05451 
.01112 
. 01508 
. 09069 
. 32523 
. 62984 
. 90347 


-0. 67876 

-. 04188 
. 01153 
. 06607 
. 23591 
. 32278 
. 14340 
. 00829 
. 04812 
. 16811 
. 24758 
. 12665 


-0.37499 

. 04980 
. 00482 
. 01494 
. 08995 
. 32481 
. 62627 
. 84821 


-0.65241 

-.04055 
.01117 
. 06344 
. 22943 
. 32360 
. 15885 
. 00801 
. 04643 
. 16253 
. 25009 
. 13898 


-0. 38284 

. 04621 
-.00029 
. 01478 
. 08913 
. 32275 
. 62148 
. 80637 


-0. 62878 

-.03938 
. 01085 
. 06119 
. 22333 
. 32468 
. 17377 
. 00775 
. 04492 
. 15750 
. 25088 
. 15150 


-0. 38910 

.04344 
-.00451 
. 01461 
. 08823 
. 31985 
. 61579 
. 77481 


-0. 60745 

-.03831 
. 01056 
. 05928 
. 21755 
. 32582 
. 18825 
.00752 
.04353 
. 15299 
. 25039 
. 16406 


-0.39403 


Eigenvector: 


04128 


2 


— .00804 


3 


. 01445 


4 


. 08725 


5 


. 31658 


6 


. 60939 


7 


. 75126 


8 




9„ 












10 













11 
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Table 3. 


Summary of results 






Reia.u.) 


^«(a.u.) 




Calculated 


Experimental 




BHr 


2.35 
2.05 
1.94 
2.01 
1.73 


a 2. 37 
b 2. 062 
cl.95 
dl.95 


-26.946 


CH4. 


—40. 167 


NH4+ 


-56.504 


CH2, 3S-g 

H3 


-38.889 
-1.593 










R 


Be 






2.1 
1.95 
1.8 
1.9 


Calculated 


Experimental 




CH2, lAi 

NH2 


105°27' 
107°29' 
106°56' 
119°23' 


dl03° 

e 103°23' 
f 104°27' 


-38. 803 
-55.454 


H2O 


—75.845 


H2O+ - 


-75.447 













a P. T.'Ford and R. E. Richards, Disc. Faraday Soc. 19, 230 (1955). 

b D. P. Stevenson and J. A. Ibers, J. Chem. Phys. 33, 762 (1960). 

c H. S. Gutowski, G. E. Pake, and R. Bersohn, J. Chem. Phys. 
2^,643 (1954). 

d G. Herzberg, Proc. Roy. Soc. (London) A262, 291 (1961). 

e K. Dressier and D. Ramsay, Phil. Trans. Roy. Soc. (London) 
A251, 552 (1959). 

f G. Herzberg, Infrared and Raman Spectra (D. Van Nostrand 
Co., New^York, 1945), p. 488. 



Table 4. Deduction of hound to correlation energy 
(a.u.) 





BH4- 


CHz, 3S-, 


H3 


Calculated Hartree-Fock 
approximation 


-26. 946 

a -.550 
a -.034 
-26.680 
e -.005 

-26. 529 

-.417 

-.319 


-38. 895 

b-.283 

-.02 

-38.855 

«-.012 

-38.689 

-.200 

-.28 


-1.593 


Heat of atomization 

Zero point energy 


C-.154, f-.151 
d -.009 


Energy of separated atoms.. 
Relativistic correction 

Hartree-Fock energy for 
separated atoms _ _ _ 


-1.5 
— 1.499 


Hartree-Fock heat of 
atomization+ zero-point 
energy., 


-.094 


Upper limit to molecular 
correlation energy 


-.070, -.067 



a A. P. AltshuUer, .T. Am. Chem. Soc. 77, 5455 (1955). 
b H. Prophet, J. Chem. Phys. 38, 2345 (1963). 
c A. Farkas and L. Farkas, Proc. Roy. Soc. A152, 124 (1935). 
d 1. Shavitt, J. Chem. Phys. 31, 1359 (1959). 

e C. W. Scherr, J. N. Silverman, and F. A. Matsen, Phys. Rev. 
128,2675 (1962). 
i R. E. Weston, J. Chem. Phys. 31, 892 (1959). 



Table 5. Partition of molecular energy, -ei 



Orbital 


la, 


2ai 


3ai 


I62 


I61 


H20,i2=1.7 












105° 


53. 6388 


9. 3362 


7. 5375 


7. 4739 


7.6305 


115° 


53. 6302 


9. 3115 


7. 5286 


7. 5075 


7.6183 


125° 


53. 6209 


9. 2866 


7. 5233 


7. 5389 


7. 6050 


135° 


53.6111 


9. 2623 


7.5231 


7. 5669 


7. 5910 


H2O, 22=1.8 












95° 


53. 5918 


9. 2388 


7. 4249 


7. 2463 


7. 6022 


105° 


53. 5838 


9. 2185 


7. 4223 


7. 2693 


7. 5904 


115° 


53. 5742 


9. 1965 


7. 4209 


7. 2940 


7. 5767 


125° 


53. 5638 


9. 1734 


7. 4241 


7. 3178 


7. 5618 


135° 


53. 5525 


9. 1501 


7. 4329 


7. 3391 


7. 5455 


150° 


53. 5351 


9. 1162 


7. 4559 


7. 3642 


7. 5205 


H2O, R=2 












105° 


53. 4903 


9. 0432 


7. 1825 


6.9110 


7. 5232 


115°._,.. 


53. 4797 


9. 0243 


7. 1972 


6. 9198 


7. 5080 


125° 


53. 4674 


9.0030 


7. 2169 


6. 9295 


7. 4905 


135° 


53. 4531 


8. 9792 


7. 2449 


6. 9373 


7. 4702 


150° 


53. 4299 


8. 9413 


7.3034 


6. 9435 


7. 4369 


NH2,JR=1.95 












95° 


41. 0537 


7. 2882 


5. 8303 


5. 7860 


6. 0372 


105° 


41. 0469 


7. 2651 


5.8117 


5. 8195 


6. 0264 


115° 


41. 0395 


7. 2420 


5. 7955 


5.8515 


6. 0146 


CH2,i?=2.1 












90° 


30. 1910 


5. 6279 


4. 4646 


4. 5631 




100° 


30. 1860 


5. 5980 


4. 4390 


4.6038 




105° 


30. 1836 


5. 5837 


4. 4267 


4. 6229 




110° 


30. 1810 


5. 570O 


4. 4148 


4. 6412 





3. Discussion 

It is evident that the approximate Gaussian calculations are adequate for the determina- 
tion of the molecular geometries of first-row hydrides. The greatest discrepancy between the 
calculated and experimental equilibrium angles is about 5° and the largest error in the equilib- 
rium internuclear separation is 0.05 a.u. This accuracy is of the order of that reported in 
reference 5, although a more accurate SCF result for H2O presented by Moccia yields slightly 
better agreement between the calculated and experimental angle. 

The success of the numerical experimentation permits confidence to be placed in the pre- 
dicted angle for H2O+, a rather large 119°. This value is totally unexpected from the pre- 
sumed nonbonding character of the Ibi orbital. Experimental evidence is ambiguous. The 
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first Rydberg excited state of H2O has been analyzed (18). Only a small change in angle was 
observed but it was accompanied by an appreciable increase in the 0-H distance. It is usually 
assumed that the first member of an up Rydberg series does not penetrate the core and iiifliHMice 
the bonding, but this is not certain. 

From Koopman's theorem the vertical ionization potential for H2O, 0.465 a.u., is approxi- 
mated by the orbital energy 0.500 a.u. However, as noted by Clementi (19), tlie ionization 
potential is given by 

where E'^^ is the reorganization energy in computing the SCF enci-gy of the ion, {E^—E^) is 
the difference in correlation energies between the neutral and ionic species, (Er—E^) is the 
difference in relativistic energies between the neutral and ionic species, and /invert is tlie int(M'nal 
energy of the ion attendant upon a vertical transition. 

(Er—E'^) is very close to zero and for our purposes we sliall assume that EVert is also snndl 
enough to be neglected. Then the change in correlation energy can be determined from Er^j 
which is approximately —0.100 a.u., to be about —0.065 a.u. This result is in good agree- 
ment with dementi's [20] esthnate of —0.074 a.u. for the correlation energy of a pair. In fact, 
if we use a better esthnate of the CH4 correlation energy [21], the correlation energy per bond- 
ing pair in CH4 is about —0.067 a.u. or less. TJie fortuitous cancellation of the reorganization 
energy and the change in correlation energy is the basis for the success of the application of 
Koopman's theorem to the prediction of ionization potentials. 

The results for H2O were obtained with SCF calculations some 2 eV higher than those re- 
ported by Moccia and some 6.5 eV from the Hartree-Fock limit estimated by Moccia. Agree- 
ment of estimated and computed pair correlation energies must be considered in tliis liglit. 

The results for H3 are the most accurate of all those reported. The asymptotic error for* 
three H atoms computed with a five iS-type Gaussian basis is 0.02 eV and the error in tlie mole- 
cule is probably less than 0.05 eV. Addition of another 7>>-type Gaussian improved the total 
molecular energy by less than 0.01 eV. 

Unfortunately, the Hartree-Fock H3 surface is of little value for understanding the kinetics 
of the H+H2 system since a radical reaction is primarily a correlation problem. There are 
also complications due to the inaccurate asymptotic behavior of the molecular orbitals. An 
estimate of the correlation energy cannot be made from the isoelectronic atoms; the intershell 
correlation is much larger in H3 than in Li due to the appreciable differential overlap in the 
molecule. 

The results for BH4~ are analogous to those reported for CH4 and NH4+ and little can be 
added to the discussion of BH4~. It should be noted, though, that the experimental heat of 
atomization is not accurately known and there is a corresponding uncertainty in the bound for 
the inolecular correlation energy. 

Evidently, the accuracy for the ^Sj state of CH2 is not comparable to that for the tetra- 
hedral molecules or H3. This is evidenced by the large correlation energy limit in table 4. 
Although the heat of atomization reported for CH2 may still be too high, the calculated inter- 
nuclear separation indicates the SCF result for CH2 is still appreciably above the Hartree-Fock 
limit. For strongly bonded molecules, it would be expected that most excited configurations 
required to represent the exact function would have larger equilibrium internuclear separa- 
tions. The Hartree-Fock equilibrium internuclear separation, therefore, should be snndler 
than the experimental value. The inadequacy of the Gaussian basis for CH2 can be attributed 
to two main defects: (1) no variation of the orbital exponents was attempted, and (2) 7>>-type 
functions are required on the H atoms. 

The final connnent is in many ways the most important. A comparison of table 5 with the 
similar data presented by Coulson and Neilson shows how ill-founded was their expectation 
that improved molecular wavefunctions would not change their results qualitatively. The 
behavior of e^ai i^^ the neighborhood of 105° is not at all shnilar to that exhibited by the curve 
for 3ai in figure 3 of reference 3. There is little point in examining the differences in detail; 
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we only need to note that the partitioned energies d reflect the influence of the other orbitals 
in such a way as to render a simple interpretation of these quantities very difficult. The most 
striking evidence of this effect is the large variation in the eiai and 61 ^i even though these molecu- 
lar orbitals vary only slightly with angle. There is little support for a one-electron type par- 
tition of the total energy from the present results. 



I thank Dr. S. Huzinaga for furnishing me with details of the accurate atom results he 
has obtained with Gaussian bases. I am also grateful to the Computer Science Center of the 
University of Maryland for permitting the use of their facilities. 
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